Heat production and energy balance in nanoengines driven by time-dependent fields. 
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We present a formalism to study the heat transport and the power developed by the local driving 
fields on a quantum system coupled to macroscopic reservoirs. We show that, quite generally, two 
important mechanisms can take place: (i) directed heat transport between reservoirs induced by the 
ac potentials and (ii) at slow driving, two oscillating out of phase forces perform work against each 
other, while the energy dissipated into the reservoirs is negligible. 
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PACS numbers: 72.10.-d,73.23.-b,73.63.-b 



The understanding of the heat transport at the micro- 
scopic realm has attracted the attention of theoreticians 
for several years now. Several studies investigate this is- 
sue in the framework of one-dimensional lattice models 
of interacting classical oscillators. [H, d, 0, 0, HI Nowa- 
days, the technological trend towards the fabrication of 
nanosize electronic devices, is boosting the theoretical 
interest in quantum transport in a variety of setups and 
materials. Recently, there have been efforts to address 
the related problems of energy transport and heat dissi- 
pation in these small-size systems, [a, 0, HI, & E3, Ell E3 

Electronic quantum transport through mesoscopic sys- 
tems has been traditionally analyzed as a response to dc- 
voltages. There are, however, alternative possibilities to 
induce net transport by using time-dependent fields as 
the generating source. Interesting examples of this kind 
have been recently realized experimentally. [lU, [lij An 
important characteristic of these systems is that directed 
motion is realized by pure ac forces thanks to the conve- 
nient breaking of relevant symmetries. [HI, [l6|, [ijj 

Energy transport in stationary conditions is achieved 
as a response to temperature and/or chemical potential 
gradients. The application of time dependent fields can 
induce net particle transport between reservoirs at the 
same chemical potentials, while it brings about by itself 
heating of the sample. Then, it is possible that ac forces 
can also generate directed heat transport between those 
reservoirs even if they are at the same temperatures. If 
so: What determines the direction of that heat flow? In 
general: What is the detailed energy balance for such a 
system? Is it possible to transport part of this energy to 
develop work? The theoretical study of the underlying 
physical processes demands a full quantum-mechanical 
treatment of the problem to evaluate the heat currents 
through the different parts of the device, as well as the 
calculation of the powers developed by the external fields. 
The aim of this work is to present a theoretical approach 
based in non-equilibrium Green's function that will al- 
low us to address details of the energy balance in the 
framework of an exactly solvable model of an electronic 
quantum pump. [18] 
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FIG. 1: (Color online) Sketch of the device. The sites of the 
chain described by H sys (t) are represented by small circles. 
The local harmonically time-dependent fields Vj(t) develop a 
power Pj. The right box encloses a unit of chain described 
by H l ' l+1 (t) where energy conservation is considered, while 
the left one corresponds to a unit described by ff( cont \ The 

— E 

arrows indicate the dc energy currents J t flowing into a bond 
(I, I + 1) of the system, J ta entering the contacts and J a 
entering the reservoirs. 



The device can be described in terms of a Hamilto- 
nian with three pieces representing the electronic driven 
system, the contact to the reservoirs and the reservoirs: 
H(t) = H s y s (t) + H cont + H res . For simplicity, we con- 
sider a two terminal setup with left (L) and right (R) 
reservoirs, and a AT-site one-dimensional lattice model 
for the driven system. We assume that the latter has 
an energy profile £° and nearest neighbor hopping el- 
ements wi. At M lattice positions, ac potentials of 
the form Vj(t) = Vj cos(Oo* + fj) ar e locally applied. 
This system Hamiltonian can, thus, be expressed as 
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H sys^ = Y,^ 1 HV- l+1 \t), being 

Jl(i.i+i)( t ) = e ,(t) c tQ +ei+1 (t) c t +lCm 



with e,(t) = 0(e? + Ejli bein S = 1/2 for 

1 < I < N and £j a = 1, for the sites that intervene in the 
connection to the reservoirs l a = 1,JV. The contacts are 
represented by hopping terms between the reservoirs and 
the latter positions: H cont = J2 a k a w a( c l a c i a + H.c). 
The reservoirs are described by free-electron models: 
jjres _ ek a c\, Ck a , which are assumed to remain 
at equilibrium with well defined chemical potentials p a 
and temperatures T a , even after being connected to the 
driven structure. 

In order to define currents of energy, we can follow 
the quantum mechanical counterpart of the procedure 
carried out in Refs. 0, H| for a classical model. We for- 
mulate the equation for the conservation of the energy 
E l)l+1 (t) = stored in a bond (1,1+ 1). The 

variation in time of this mean value is calculated by re- 
course to Ehrenfest's theorem which casts: 



dE u+l {t) 
dt 

dE a ,i a {t) 



= J) 



f +1 (t)-J?(t) + l(Pl(t) + Pi + i(t)}?) 

(3) 



where the first equation corresponds to bonds of the sys- 
tem without sites in contact to reservoirs (i.e. 1 < I < 
N — 1), as the one enclosed by the right box of Fig. [TJ 
while the second one corresponds to the left bond that 
establishes the contact between system and reservoir (see 
left box of Fig. []}. The first two terms of ([2]) represent a 
discretized version of the divergence of the energy current 
flowing through the bond, while the last two terms repre- 
sent the power developed by the external forces and are 
equal to (dH (l ' l+1, >(t)/dt). In Eq. ©, the first current 
represents the flow of energy towards the system, the sec- 
ond one is the flow of energy entering the reservoir, and 
the third one, the power developed by the external forces. 

Denoting Pj.j>(t) — (c^,(t)cj(t)), the explicit expres- 
sions for the different energy currents read: 



Jl E (t) 
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\ei{t) pi t i-i(t)wi-i - pi t i +1 (t)wi (4) 
+wi-ipi-i t i+i(t)wi\, 
-2Im| w a pi ai k a (t)E ka X, (5) 



and a similar expression for J l ; (t), while the power de- 
veloped by the external forces is: 



M 



dt 



(6) 



In order to analyze the energy balance we focus on the dc 
components of the energy currents and powers done by 
the external forces. The conservation of the energy (|2|) 
and © implies dE l<l+1 (t) / dt = and dE aJa (t)/dt = 0, 
which defines continuity equations for the dc energy cur- 
rents J i , J a and powers Pi, where A = 1/tq /J" dtA(t), 
being tq = 2tt/D,q the period of the oscillating time- 
dependent fields. 

To evaluate the different energy currents and pow- 
ers we employ the treatment based in Keldysh non- 
equilibrium Green's functions of Ref. [l8[, which has 
been useful to study charge transport in the kind of sys- 
tems we are considering. The mean values of observ- 
ablcs entering the corresponding expressions can be ex- 
pressed in terms of lesser Green's functions as follows 
Pj.j'(.t) = —iGfji(t,t). The latter is evaluated from a 
Dyson equation, which for j, j' lying on the central sys- 
tem reads: 



G<(t,t) 
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a—L,R kk'— — c 



-ikflot 



du! 
2^ 



xg ltla (k + k',uj)g;, l (k',uj), 



(7) 



being T<(lu) — if a (u)T a (uj), where the Fermi func- 
tion fa(oj) depends on p a and T a , and T a (uj) = 
2nJ2k \ e k a i a \ 2 S(uj — Sk a )- The Green's function 
Gi,i'{k,oj) is the fc-th Fourier coefficient of the Fourier 
transform of the retarded Green's function, which can 
be exactly evaluated with convenient methods. [18| The 
above expression ||7J) can be used in ((4]) and ([6]) to evalu- 
ate J ; and Pu . Using properties of the Green's function 
the dc component of the energy current flowing into 
the reservoir reads: 



7 « = E E / f% + /A>)r a (w + *A))i>H 

X \Q lath (k, L0)\ 2 \f p (u) -fa(w + fcf> )l . (8) 



We are interested in the case of reservoirs at the same 
chemical potentials fi a — /i,Va. Following the same line 
as in stationary transport at linear response, ll| we define 

the de-heat current as J 1 — J 7 — p,J 1 , the latter term 
being the convective flow which depends on the dc-charge 



current J ' .[20( As this current is conserved, J ' obeys the 

— E 

same continuity equation satisfied by J 1 . We now use 
the above theoretical framework to analyze two generic 
mechanisms that can take place in quantum pumps. The 
first one concerns the possibility of achieving directed heat 
transport between reservoirs. To this end, notice, that 
the dc-heat current along the lead a can be splitted as 
the addition of a generated and a pumped contribution: 

= J^} 3 ^ + J^} P \ where the pumped component reads 
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FIG. 2: (Color online) Pumped heat current (left) and total 
power P = Pi + P2 developed by the forces (right) as func- 
tions of temperature and the pumping frequency. Parameters 
are N = 22, E = 1, V = 0.2, T = 0, tp = ty/2 and /i = -1.5. 



FIG. 3: (Color online) Total heat current as functions of T 
for different pumping frequencies. Thin and thick lines corre- 
spond to and J%, respectively. Parameters are as in Fig. 
[21 The values of fio are indicated in the plots. 



explicitly: 
-Q( P ) f du) 



J J( W - M ){2r<Hlm[^, iQ (0^)]+ (9) 

OO 

E E T a {uJ + k£l )T<{u)\g^ h {k,u)\ 2 ). 



/3=L,Rk=-oo 



This component can be proved to satisfy 

Sa=L r J < c} P = 0, meaning that heat can be ex- 
tracted from a given reservoir and injected into the 
other one. From the continuity equation for the dc 

heat current, it follows that J2f=i Pi 3 = J2 a =L.R^a 9 j 

indicating that the contribution accounts for the 

heat generated by the external forces, which is dissipated 
into the reservoirs. 

An example of the behavior of these two different con- 
tributions to the heat flows at the reservoirs is shown in 
Fig. [2] for a two barrier setup in contact to reservoirs 
at the same temperature Tl — Tr — T. We consider 
£;° = E, lj = 2, N - 1 and ef. = 0,1 j= lj, with two 
oscillating potentials with the same amplitude V and a 
phase-lag tpi = tp,<^2 = applied at the barriers. Un- 

Q 

der such conditions a dc charge current J is induced 

which behaves like J oc V 2 sm(tp) at small V. [IH The 
pumped component of the heat current, shown in the 
left panel of Fig. [2l flows outwards in a reservoir and in- 
wards in the other one. It exhibits a complex structure of 
maxima, minima and sign inversions as a function of CIq. 
The details of these features are model-dependent. They 
are consequence of the electronic propagation through 
a structure with discrete energy levels and quantum in- 
terference that takes place when the pumping frequency 
£7o is resonant with the energy difference between these 



levels. Q The generated heat currents at the differ- 
ent leads J L R , also display a complex landscape in the 
T, Qq plane. Their sum is, however, always positive and 
equals the total developed power P = P\ + P 2 , which 
is a monotonous increasing function of il and decreases 
with T as shown in the right panel of Fig. [2l At T = 
dissipation dominates and masks the heat pumping ef- 
fect. For our example, the total heat currents J L , as 
functions of T are plotted in Fig. [3] There it is seen that 
the sign of both currents is positive at T = 0, which indi- 
cates that the flow goes from the central system towards 
the reservoirs. Such a behavior is the one expected from 
considerations based on general thermodynamics since, 
at T = there is no heat at the reservoirs amenable to 
be transported. However, the results of Fig. [3] show that 
at finite T a regime exists for low pumping frequencies 
Slo where heat pumping takes place. In fact, the signs of 
the heat currents are different, indicating that they leave 
one of the reservoirs and enter the other one. For higher 
f2o dissipation is again dominant and heat flows always 
from the central system into the two reservoirs. 

The second remarkable mechanism we would like to 
analyze is the possibility of extracting work from the sys- 
tem. This issue can be addressed on the basis of a per- 
turbative solution of the Dyson equation which allows for 
the evaluation of P{. at the lowest order in Vj . [13, [3 It 
casts: 
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FIG. 4: (Color online) Main Panel: Heat currents J l , Jr 
(dashed thin and thick lines) and Jg (circles) along the left 
and right leads and the central system, respectively for^t = 
— 1.5. Powers developed by the external fields P\,P2 (thin 
and thick solid lines). Inset: Pi, P2 as function of /j for Oo = 
0.01. Parameters are as in Fig. [5] 



where the coefficients are 



A 



(i) 



X (2) 



-OO 
+ OO 



^^{7^(^)7^.(0;)}, 

^Bfi{ 7i i(a;)7^H}, 
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(11) 



with 7j ,m = Ea=£,H^H[G? J ,i«H]*^,i«H and 

7 ± (w) = [g° ^ (w + fi ) ± G° >(i (w - Oo)] , being G°,, (w) 
the equilibrium retarded Green's function of the sys- 
tem in contact to the reservoirs but without the time- 
dependent fields. In the low frequency limit, Ay oc J7q 

(2) 

while AL oc £Iq ■ Thus, this solution indicates that quan- 
tum coherence in the wave function propagation along 
the structure, which rules the behavior of the charge cur- 



rent, also plays a role in the way in which energy is pro- 
vided and exchanged. In particular, for more than one os- 
dilating fields, the terms oc AL/ dominate at low enough 
Oo- Since these terms can be positive for some fields and 
negative for other ones, this enables a scenario where the 
total energy P — Ylj=i ^ i ls dissipated to the reservoirs 
at a ratio oc Oq, while a larger amount of energy oc Qq is 
exchanged between the different pumping centers. Such 
an effect is, in fact, observed for some parameters in the 
example of the two-barrier setup. Results are shown in 
Fig. [4] for T — with tp = ir/2. The dc powers are 

shown along with the dc heat currents J L , J R flowing 
to the reservoirs and the one flowing within the system 

between the two barriers Jg. The exchange of energy 
between the two fields is further highlighted in the inset 
where the two dc powers P\ and P2 as functions of /i are 
shown. It is also interesting to note that the direction 
of the heat flow between the two pumping centers, goes 
from the field doing the largest power towards the other 
one. Instead, in the reservoirs the heat currents are al- 
ways positive, indicating that they flow inwards. These 
features are in line with the idea that the fields heat lo- 
cally the sample inhomogeneously, then the heat current 
flows from the hottest to the coldest regions. 

To conclude, we introduced a treatment based on the 
local conservation of the energy in order to investigate 
details of the energy transport in open quantum systems 
driven by time-dependent fields. We identified two in- 
teresting mechanisms like the possibility of achieving di- 
rected transport of heat between reservoirs at finite but 
identical temperature as well as extracting useful energy 
to make work against the external forces. We illustrated 
these effects in an exactly solvable model of an electronic 
quantum pump. Our results for reservoirs at T = also 
suggest that the direction of the heat flows seem to be 
ruled by an effective heat gradient induced by the appli- 
cation of the external driving. 
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